Prom BCS to BEC Superconductivity: Spectroscopic Consequences 
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The evolution from BCS to BEC superconductivity in the s-wave and d-wave channels is analyzed 
at zero temperature for a two-dimensional superconductor. Spectroscopic quantities for s-wave and 
d-wave systems present fundamental differences when particle density and attraction strength are 
varied. A detailed analysis of single quasiparticle properties (excitation spectrum, momentum dis- 
tribution, spectral function and density of states) indicates that the evolution of these spectroscopic 
quantities in the d-wave case is not smooth, unlike the situation encountered for the s-wave system. 
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The problem of the evolution from BCS to BEC su- 
perconductivity is an old one |l],|| but recently it has 
received considerable attention in connection with high 
temperature superconductors, f|-|j| where strong devi- 
ations from the BCS behavior have been observed ex- 
perimentally in spectroscopic quantities |jTo| at low tem- 
peratures. Furthermore, the recent discovery of BEC in 
atomic systems and the exciting possibility of BEC in 
atomic Fermi systems raises the question of the evolution 
from BCS to BEC in atomic Fermi systems as well. |L2| 

In this work we address the question of whether the 
evolution of spectroscopic quantities from a BCS to a 
BEC superconductor is smooth at zero temperature. For 
this purpose we study the single quasiparticle properties 
(excitation spectrum, momentum distribution, spectral 
function and density of states) as a function of attraction 
strength or particle density for the s-wave and d-wave 
cases. In anticipation of the main results, we must say 
that the evolution of spectroscopic properties in the s- 
wave case is smooth, while in the d-wave it is not. The 
main reasoning for this statement is as follows. 

Quite generally the evolution from BCS to BEC su- 
perconductivity can be characterized by two parame- 
ters: the chemical potential /i and the Cooper pair size 
£, P air- The BCS limit is characterized by a positive chem- 
ical potential fj, = ep and a large size of Cooper pairs 
impair ^ kp ), while the BEC regime is characterized 
by a large and negative chemical potential /j, = —E- ( 



where Ei CJ is the binding energy of the two-body prob- 
lem in the £ angular momentum channel, and by a 
small size of pairs (£, pa ir kp 1 ). Here £ = (or s) indi- 
cates the s-wave channel, while £ = 2 (or d) indicates the 
d-wave channel. The excitation spectrum at zero tem- 

1/2 

perature has the form E e (k) = ((e k — fi) 2 + | A e (k)| 2 ) , 
where £k = k 2 /2m and A^(k) = Aoghg(k) cos(£<j>), with 
k = |k| . In the s-wave case the excitation spectrum E s (k) 
is gapped for all k, and it increases smoothly from the 
BCS to the BEC limit. As a result the quantities that 
depend directly on the excitation spectrum E s (k.) also 
evolve smoothly. For instance, the quasiparticle density 
of states N s (oj) at low frequencies is always zero, since 



there are no available states inside the gap. Thus, contri- 
butions from single quasiparticle excitations to thermo- 
dynamic quantities are always exponentially small at low 
temperatures. In the d-wave case the situation is qualita- 
tively different. For fj, > the superconductor is gapless 
at the Dirac points k — k^ = ^/2m/i, <fi = ±7r/4, ±37r/4, 
while for /x < the superconductor acquires a finite gap. 
The line /j, = separates two regimes with qualitatively 
different behavior. This has important consequences for 
the momentum distribution, spectral function, and den- 
sity of states. The quasiparticle density of states Nd(w) 
changes discontinously at low frequencies from linear in 
u) for /i > (where -E^k) is linear in momentum close to 
the Dirac points), to a constant at \i = (where i?d(k) is 
quadratic for small momenta), to zero for fi < (where 
E d (k) ~ \fj,\ +0(k 2 ) for small k). Thus, contributions 
from single quasiparticle excitations to thermodynamic 
quantities at low temperatures also exhibit singular be- 
havior in the vicinity of /i = 0. 

In order to analyze how the spectroscopic quantities 
change from the BCS to BEC limit, we start with the 
two dimensional Hamiltonian 



H = 



kk'q 



lL q & k', 



(1) 



where 6kq = V ; -k+q/2iV'k+q/2T- The interaction po- 
tential Vk k ' is expanded in its angular momentum 
components as Vkk' = Sfc^oo ^ki ^ ex P(^0fcfc')j where 
acos(k ■ k') is the angle betwen the vectors k 



»kk' 



and k' and V$ = 2tt / °° drrJi(kr)Je(k'r)V(r). The 
index £ labels angular momentum states in two 
spatial dimensions, with £ = 0,±1,±2, corre- 
sponding to s,p,d,..., channels respectively. A pos- 
sible choice of the real space potential can be 
V{r) = ViQ(Ri - r) - V B(r - i?i)6(i? - r), which is 
repulsive at short distances r < R\, attractive at inter- 
mediate distances R\ < r < Rq, and vanishes for r > Ro. 
Quite generally any short ranged real space potential 
V[r) with range Rq leads to a vSS which is sepa- 
rable for small momenta, provided that kRo «1 or 
k'Ro <C 1. In the simpler limit when both kRo <C 1 and 
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k'R Q < 1, V$ ~ k e k' e (2n/2 2t ) J °° drr 2e+1 V(r), where I 
is assumed to be positive for definiteness. Notice here 
that oc k e k' e , thus for the s-wave case cx const., 
while for the d-wave case v£$ oc k 2 k' 2 . In the opposite 

limit, kR > 1 or k' R Q ^> 1, the potential V^} is cer- 
tainly not separable. In the simpler limit when both 
kR ^> 1 and k' Rq 1, V^, mixes different k and k' 
and shows an oscillatory behavior which is dependent on 
the exact form of V(r), with a decaying envelope propor- 
tional to k- l l 2 k'~ 1/2 . 

Under these circumstances, quite generally it is not 
possible to find a separable potential in momentum 
space Vkk' = — Au>*(k)u>(k'), nevertheless in the spirit 
of ref. H we choose to study a separable poten- 
tial that contains most of the general features de- 
scribed above. In addition, we consider only sin- 
glet superconductivity, where the s-wave and the d- 
wave channels are studied separately. For this purpose, 
we use the separable potential Vkk' = — A^u^(k)w;^(k'). 
The interaction term we(k.) can be written as the 
product of two functions, u^(k) = he(k)gi(k), where 
h t (k) = (fc/fci)7[l + (fc/fc )] £+1/2 controls the range of 
the interaction and <j^(k) = cos(^) is the angular depen- 
dence of the interaction. Here fco ~ Rq 1 and k\ sets the 
scale at low momenta. At zero temperature, we assume 
that pairing occurs with the same total momentum q = 
only. This simplifying feature leads to the following sad- 
dle point and number equations, 

i_ = v Mk)| 2 

A, ^ 2S/(k) ' 1 1 

n = 2£)n/(k), (3) 

k 

where n^(k) = [1 — (ek — aO/^O 4 )] /2 is the momentum 
distribution, E e (k) = ((e k - fi) 2 + \A e (k)\ 2 ) 1/2 , is the 
single particle excitation energy, and (k) = Ao£Wi(k). 
is the order paremeter. For a given interaction range 
R ~ fc^ 1 , the transition from the BCS limit (largely 
overlaping pairs) to the BEC limit of (weakly overlap- 
ing pairs) may occur either by changing the attraction 
strength Xi or the density n. In either case, this evolution 
can be safely analyzed with the approximations used here 
provided that the system is dilute enough, i.e., n <C k^. 
This means that below a maximum density n max ~ 
kFm ax , the interaction range Rq is much smaller than 
the interparticle spacing kp^ ax , Rq <C kpmax' or ec l u i v ~ 
alently ko/kF max 1. Thus we choose to scale all en- 
ergies with respect to the maximal Fermi energy e^ maa ., 
which fixes the maximum density n = n rnax = 2pep max , 
and all momenta with respect to kFmax — V^ me F m ax- 
The coupling constant is scaled with respect to the two- 
dimensional density of states p. From now on we use this 



scaling. The numerical solutions for Ao^ and p, when 
k\ = fco — 10 are shown in Fig. |l|, for fixed density 
n = 1, and changing Xg. Similar plots can also be made 
for fixed interaction and varying density n. 
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FIG. 1. (a) The order parameter Ao and (b) the chemical 
potential p as a function of coupling at fixed density n = 1 
and fci = fco = 10 for both s- and d-wave channels. 

In the BCS limit the amplitude of the order parameter 
(4> = 0) is given by 

A £ (fc M ) ~ exp [2(A-/(fc M ) - Xj^/ht 2 ^)) . 

With our choice of he(k), A 0c i(fc AI ) ~ 8 + /i/24ei + 
C([/x/ei] 2 ), valid for fi/ei <C 1, where e\ = k 2 . The 
ratios between A^(fc„), and the critical temperature 
T c i satisfy the usual relations A s (fc M )/T cs = 1.76, and 
A d (fc^)/T C(J = 2.14. In Fig. |l| A Qd and // have a second 
order discontinuity as a function of A^. This discontinu- 
ity occurs when p = in both Aod and p, for varying 
interaction A^ or varying density n. The line p = for 
a d-wave system is very special as it will be seen in the 
following discussion of spectroscopic quantities. 

The first spectroscopic quantity to be analyzed is the 
single quasiparticle excitation spectrum Ee(k). Let us 
discuss first the s-wave case in the zero range interaction 
limit fco — > oo. For p > the excitation spectrum has an 
isotropic gap at fc = fc„, E g (k^) — |A s (fc„)|. This gap 
is completely isotropic in the vicinity of fc M . At the in- 
termediate regime, when p = 0, the gap takes the value 
E g (0) = |A S (Q)|, when the chemical potential becomes 
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negative towards the BEC limit, the minimum of the en- 
ergy gap remains at k = 0, E g (0) = ((j, 2 + | A s (0)| 2 ) 1/2 . 
When kg is finite the position of the minimum gap 
changes, but the excitation spectrum is always gapped. 

In the d-wave case the situation is qualitatively differ- 
ent. For fi > 0, including the BCS limit, the excitation 
spectrum is gapless at along the special directions 
<t> = ±7r/4, ±37r/4, near which the excitation spectrum 
disperses linearly with momentum. The energy gap at 
k = k ll and <f> — 0, E g {k ll ) — \Ad{k^)\ is a nonmonotonic 
function of k^ for fixed density, and thus a nonmono- 
tonic function of A^. The maximum E g (k^) is reached 
at intermediate values of fj, > 0. At fi — 0, the mini- 
mal gap is E g (0) = |Ac;(0)| = 0, and occurs at the sin- 
gle point k = 0. In this case the excitation spectrum 

1/2 

is E d (k) = (ef. + |A d (k)| 2 ) , which behaves quadrat- 
ically for small momenta at any given angle 4>, since 
A d (k) ~ k 2 cos(2(/>) and £k = k 2 /2m. The shrinking of 
the energy gap to zero at k = is a consequence of the di- 
minishing pairing interaction hd^k^) for fi — > 0. As soon 
as (i < 0, including the BEC limit, a full gap in the ex- 
citation spectrum appears, but the minimal gap remains 
at k = with value E g (0) = \fi\ since A d (0) = 0. Thus, 
the /i = line separates a gapless d-wave superconduc- 
tor (/j, > 0) from a fully gapped d-wave superconductor 
(/i < 0). Fig. § shows the lines where /j = 0on the graph 
of n vs. X(. Notice in Fig. || that the low density limit 
of the s-wave system is always Bose-like, i.e., a two-body 
bound state appears at arbitrarily small A s . On the other 
hand, the d-wave system is qualitatively different: it is 
BCS-like for A^ < X c d and Bose-like for A^ > \ c d, where 
the critical coupling A c separating the two regimes is fi- 
nite, i.e., the appearance of a two-body bound state in 
the d-wave case requires finite A^- 



/j > 0. In the vicinity of k^ the momentum dis- 
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FIG. 2. The line = for both s- and d-wave order pa- 
rameters for n — 1 and fci = ko = 10. 

The behavior of the excitation spectrum as a func- 
tion of /i has important consequences on the momen- 
tum distribution ni(k) at zero temperature. In the s- 
wave case the momentum distribution is isotropic in 
momentum space. Here we discuss briefly the be- 
havior at low momenta for three different regimes: 
fx > 0, /i = 0, and fi < 0. Let us start with 



2fc M $fc/A a (fc„)]/2. At 
[l+7 P (l + afc/2fc )]/2, 
and a = A 2 /( M 2 + A 2 J. 



tribution is n s (k^ + 5k) ~ [1 
low k it behaves as n s (k) : 

where j p = /x/V V + Aq s 
When fi = 0, the momentum distribution at small 
momenta is n s (k) ~ (l — fc 2 /Ao s ) /2. For negative 
fi, n s (k) = [1 — 7n(l + ak/2ko)] /2 for small k, with 
In = l^l/vA^ + Aos- Notice that n s (0) is a continuous 
function of fi. In fact n s (k) is a smooth function of [i for 
all momenta. This is not the case for a d-wave system, 
which shall be discussed next. 




FIG. 3. The momentum distribution of quasiparticles for 
4> = 0, n = 1, fci = ko = 10, and several values of /x for a 
d-wave order parameter. The inset shows results for /i < 0. 

The momentum distribution in in the d-wave case is 
anisotropic, having the form rid(k) = [l — sgn(fc 2 — /_*)] 
along the direction of the nodes {<j> = ±7r/4, ±37r/4). 
This behavior already signals discontinuity of Ud(k) as 
a function of fi at k = 0, a suspicion further confirmed 
by analyzing the more interesting direction 4> — and its 
equivalents <f> = ±7r/2,7r. Near k^ the momentum dis- 
tribution is nd(kfi + 6k) ~ (1 — 2k fl 6k/ Adikft)) /2. On 
the other hand, the momentum distribution behaves as 
n d (fe)^l-(Ay M 2 )(fcV4^) 

for small momenta. When /i = the momen- 
tum distribution at k = is nd(0) — (1 — k) /2, where 
k = (1 + A 2 )d /kf)~ 1 / 2 . When fi < 0, then 



n d (k) ~ (A 2 d /^ 2 )(fc 4 /4fcf) for small k. Notice the dis- 
continuity of the momentum distribution at low fc, when 
chemical potential crosses zero. This discontinuity, which 
is illustrated in Fig. [| coincides with the collapse of the 
four Dirac points to a single point at fc M = 0, and with 
the appearance of a full gap as soon as fi < 0. 

The qualitative changes in -E^(k) and ne(k), as a func- 
tion of fj,, affect substantially the quasiparticle density of 
states N £ (uj) = N^ +) (uj) + N^~\u), where 

iV, (+) M = [ d 2 k[l - n e (k)} S(uj - E e (kj), (4) 



corresponds to adding a quasiparticle, and 
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N^\ui) = [27T]- 1 I d 2 kne(k)5(iu + E e {k)), 



(5) 



corresponds to removing a quasiparticle. In the s-wave 
case N s (oj) is always zero at low frequencies, since the 
excitation spectrum is gapped for all \±. On the other 
hand, Nd{u>) changes discontinously at low frequencies, 
from linear in w for /i > 0, where Ed(k) is linear in 
momentum close to the nodes, to a constant at /i = 
(where ^(k) oc fc 2 at low fc), to zero for n < (where 
Ed(k) ~ |/x| +0(k 2 ) for small fc), as can be seen in Fig. 0. 
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FIG. 4. Density of states for a d-wave order parameter near 
(i — 0, for n = 1, fci = kg = 10, and varying A^. Notice the 
strong asymmetry (absence of quasiparticle-quasihole symme- 
try) between uj > and uj < 0. 

Lastly, it is important to point out qualitative differ- 
ences in thermodynamic quantities, e.g. specific heat C 
and spin susceptibility %, during the evolution from BCS 
to BEC superconductivity at low temperatures. The 
contributions from quasiparticles to C and \ are expo- 
nentially small at low temperatures in the s-wave case 
for all fj,, given that the excitation spectrum is always 
gapped. The situation is qualitatively different in the d- 
wave case where C oc T 2 , and x « T for /i > 0; C « T, 
and x const, for fi = 0; C oc T -1 exp(— |/i|/T), and 
X oc exp(— \n\/T) for /i < 0. 

In summary we studied the low temperature evolution 
from BCS to BEC superconductivity for varying density 
and interaction strength in both s-wave and d-wave chan- 
nels. In the s-wave case the excitation spectrum is always 
gapped, and the momentum distribution is a continuous 
function of //. However, in the d-wave case the excitation 
spectrum is gapless for /i > and acquires a full gap for 
(i < 0. Furthermore, the momentum distribution is dis- 
continuous at low fc, as [i crosses zero. As a result, the 
changes in spectroscopic and thermodynamic properties 
near — are dramatic at low temperatures. Jig] The 
line y, = in the n vs. plane seems to correspond to 
a quantum critical line. 
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